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RESUMEN 

Una buena parte del Centra Norte de Quito se halla sobre el segmento de 
falla llumbisi-La Bota, ILB, donde se espera un sismo de magnitud 6.2, con un 
periodo de recurrencia que esta alrededor de los 600 anos, con foco superficial, lo 
que da una gran peligrosidad sismica. Por este motivo, se realizo un estudio 
donde se definieron cuatro microzonas y para cada una de ellas se determino 
espectros de respuesta cuyas ordenadas espectrales para la zona de penodos 
cortos son mas altas que las que se obtienen con el espectro de la Norma 
Ecuatoriana de la Construccion, NEC-11 o con el espectro hallado con los factores 
de sitio de ERN-12. 

Pero es solo mas alto en la zona de penodos cortos y mas bajo en la zona 
de penodos intermedios y largos de tal manera que al analizar una estructura con 
el espectro asociado a la falla ILB (sismo cortical) se esta sub dimensionando la 
peligrosidad sismica para estructuras de mediana y gran altura. Esto se lo ilustra 
con el analisis sismico de una estructura de acero de 6 pisos, aplicando el Metodo 
del Espectro de Capacidad. 

Por otra parte, se presenta con bastante detalle el marco teorico del 
Metodo del Espectro de Capacidad y se describe el uso del sistema de 
computacion CEINCI-LAB para encontrar: la curva de capacidad sismica 
resistente aplicando analisis estatico no lineal; el espectro de capacidad; el punto 
de desempeno y los desplazamientos y fuerzas laterales en cada uno de los pisos. 

ABSTRACT 

Much of North Central Quito is on the fault segment llumbisf-La Bota, ILB, 
where an earthquake of magnitude 6.2 is expected recurrently, approximately 
every 600 years, with shallow focus, which thus produces a great seismic hazard. 
Therefore, a study was carried out in which four microzones were defined and for 
each of them the spectral ordinates of the response spectra was determined. For 
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the area of short periods they are much higher than those obtained with the 
spectrum of Standard Ecuadorian Construction, NEC-11 or the spectrum found 
with site factors ERN-12. 

But it is only higher in the area of short periods and lower in the area of 
intermediate and long periods so that when analyzing a structure with the spectrum 
associated with fault ILB (cortical earthquake) there is a seismic danger for medium 
to high structures. This is illustrated with the seismic analysis of a six storey steel 
structure, applying the Capacity Spectrum Method. 

On the other hand, the theoretical framework of the Capacity Spectrum 
Method is presented in very much detail and describes the use of the computer 
system CEINCI-LAB to find: The seismic resistant capacity curve using nonlinear 
static analysis; spectrum capacity; the point of performance and lateral 
displacements and forces on each of the storeys. 

1 INTRODUCCION 



En la fotograffa de la figura 1, se aprecia parte del Centra Norte, en primer 
piano y al fondo el Sur de Quito, al lado derecho se aprecia la Cordillera 
Occidental y en la parte izquierda a todo lo largo se ven las lomas de Puengasf al 
Sur y las lomas de llumbisf-La Bota, ILB. 




Figura 1 Panoramica de Quito y dos segmentos de las fallas ciegas de Quito: el 
segmento de Puengasf al Sur y el llumbisf-La Bota en el Centra Norte. 
Fotograffa: Trujillo (2014). 

Estas lomas constituyen dos de los segmentos de las fallas ciegas 
inversas de Quito que se estan levantando entre 3 y 4 mm., por aho debido a las 
fuerzas de compresion que se dan en la corteza terrestre, fruto del movimiento 
tectonico de la placa de Nazca y la placa Americana. 
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El segmento de falla ILB tiene un angulo de buzamiento que esta 
alrededor de los 55° hacia el oeste, y el piano de falla tiene un ancho que varia 
entre 10 y 12 km. (Alvarado et al. 2014). La maxima magnitud que se espera en 
este segmento de falla es M w = 6.2 la misma que fue obtenida aplicando las 
ecuaciones de Leonard (2010) para una longitud de ruptura de 15 km., igual 
magnitud se encuentra al considerar un area de ruptura de 176 km 2 . Alvarado et 
al. (2014). 

Se define el penodo de recurrencia como el tiempo que transcurre entre la 
ocurrencia de dos eventos de las mismas caracteristicas. Existen varios modelos 
para hallar este penodo uno de ellos es el denominado Terremoto Caracteristico 
propuesto por Wesnousky, (1986). Con este modelo se encontro que el penodo de 
recurrencia para un sismo de magnitud 6.2, asociado al segmento ILB es de 138 
anos. Alvarado etal. (2014). 

El modelo de Terremoto Caracteristico considera una periodicidad 
perfecta, para el caso de ILB se tendria que durante 138 se va acumulando la 
energfa producto de la tasa de movimiento de la falla y luego se genera el sismo 
de magnitud maxima 6.2. Pero en el intervalo del penodo de recurrencia se 
producen sismos de menor magnitud y existe una liberacion de energfa la cual no 
es reconocida en el modelo de Wesnousky (1986). 

Un modelo que si considera la ocurrencia de sismos de menor magnitud, 
para hallar el penodo de recurrencia, es el propuesto por Gutenberg y Richter (G & 
R) pero considerando un doble truncamiento de la magnitud desde un valor 
mfnimo M min hasta un valor maximo M max (Cosentino et al. 1977). Con este 
modelo se tiene: 

/( J n) = -^«= P e-«-» (1) 

M mln 

Donde f{m) es la funcion de densidad de la magnitud; N Mmin es la tasa 
acumulada de sismos de magnitud mayor o igual a M min ; p es el parametro de la 
formulacion clasica del modelo de G & R. En Rivas et al. (2014) esta muy bien 
detallado el calculo del penodo de recurrencia de cada uno de los segmentos de 
las fallas ciegas de Quito, aplicando el modelo de G & R modificado; los resultados 
se presentan en la Tabla 1, para diferentes rangos de magnitud y para los 
segmentos de falla de Puengasf; llumbisf-La Bota ILB; Carcelen el Inca, CEI; 
Bellavista-Catequilla, BC y Tangahuilla. 



Tabla 1. Penodo de Recurrencia, hallados con modelo de Gutenberg y Richter 
Modificado. (Rivas etal. 2014) 



Rango de 
magnitudes 


Periodo de recurrencia (anos) 


PUESGASI 


ILB 


CEI 


BC 


Tangahuilla 


[5,0 - 5,5) 


20-35 


18-30 


27-39 


18-31 


23-34 


[5,5 - 6.0) 


62-87 


56-75 


85-130 


58-78 


65-97 


[6,0 < 


164-262 


179-279 




169-279 


579- 1016 


Mmax 


1224-2190 

(Mw6,4) 


610-981 

(Mw6,2) 


549 - 952 

(Mw5,9) 


908- 1630 

(Mw6,3) 


579- 1016 

(Mw6,0) 
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Para el segmento ILB, el periodo de recurrencia de sismos de magnitud 
entre 5 y 5.5 esta entre 18 y 30 anos; para sismos entre 5.5 y 6 esta entre 56 y 75 
anos; para sismos de magnitud mayor a 6 el periodo de recurrencia esta entre 179 
y 279 anos. Finalmente para el sismo de magnitud maxima esperada de 6.2, el 
periodo de recurrencia esta entre 610 y 981 anos. 

Si se considera que en los ultimos cuatro siglos no ha habido en Quito un 
sismo de magnitud 6 o mas asociado a los segmentos de fallas ciegas y 
considerando que estan en movimiento continuo, quiere decir que hay una 
cantidad de energfa acumulada muy significativa que puede liberarse en cualquier 
momento ya que se encuentra dentro de los penodos de recurrencia esperados 
(Ver tabla 1). Se destaca que el modelo de G & R modificado es conservador con 
respecto al modelo de Terremoto Caracteristico. De ahi nace la necesidad de 
conocer las formas espectrales de los sismos que se producirian en cada uno de 
los segmentos; para el Centra Norte de Quito, el segmento de falla mas proximo 
es el de ILB. 

2 ESPECTROS PARA EL CENTRO NORTE DE QUITO 

En la figura 2 se presentan las parroquias del Centra Norte de Quito, 
destacando que la parte superior de las Lomas de llumbisi-La Bota se encuentran 
al lado derecho de la: Mariscal Sucre, lhaquito y Jipijapa; estas y las restantes que 
se hallan a la izquierda se encuentran en el bloque levantando (Hanging Wall), en 
cambio el sector de Cumbaya y los valles se encuentran en el bloque deprimido 
(Foot Wall). 

En la figura 2 se presenta la malla de puntos en los cuales se hallo el 
espectro que se generaria debido a un sismo de magnitud 6.2, con hipocentro 
ubicado en el centra de gravedad del piano de falla ILB. Los espectros fueron 
obtenidos utilizando los modelos de movimientos fuertes desarrollados por 
Campbell y Borzognia (2013); Abrahamson, Silva y kamai (2013); y el modelo de 
Zhao etal. (2006). 

En Aguiar et al. (2014) estan muy bien detallados la base de datos 
(acelerogramas utilizados), ecuaciones de calculo con las tablas de los 
coeficientes para diferentes penodos y el alcance de cada uno de los modelos. Sin 
embargo de ello se indica que los modelos C & B (2013) y el ASK (2013) permiten 
encontrar espectros de respuesta para 5% de amortiguamiento asociados a fallas 
corticales (caso de las fallas ciegas de Quito). En cambio el modelo de Zhao et al. 
(2006) permite hallar espectros asociados a tres fuentes sfsmicas, a saber: 
cortical, interface e intraplaca. 

Una vez que se encontro los espectros en cada punto de la malla y con 
cada uno de los modelos se obtuvo un espectro ponderado de la siguiente manera 
35% del espectro de C & B (2013) mas el 35% del espectro de ASK (2013) mas el 
30% del espectro de Zhao etal. (2006). 

Se pondero de la forma indicada debido a que el valor de la incertidumbre 
que se halla con el modelo de Zhao etal. (2006) es mayor al que se encuentra con 
los otros modelos. 



5 Peligrosidad sismica de Quito y el Metodo del Espectro de Capacidad 




Figura 2 Centra Norte de Quito y malla donde se obtuvieron espectros de 
respuesta elastica asociados a un sismo de magnitud 6.2 y foco superficial en el 
segmento de falla ILB. (Aguiar etal. 2015) 




Piano de ruptura (Mw=6.2) 
del segmento llumbisi'-La Bota 

• Sitios de estudio 



ADtfllNJ 

Parroquias del Centra Norte de Quito 
Administracion Centra Norte 
Administraciones zonales 



Zona 1: valores maximos HW 
I I Zona 2: valores moderados HW 
I I Zona 3: valores moderados FW 
■ Zona 4: valores mi'nimos HW 



Figura 3 Microzonificacion del Centra Norte de Quito 
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Con los espectros hallados se obtuvieron las cuatro microzonas del Centra 
Norte de Quito, indicadas en la figura 3 y para cada una de ellas se encontro el 
espectro de mayores ordenadas espectrales para un 84% de confiabilidad. Antes 
de presentar los espectros se destaca que la parte achurada de la figura 3 
corresponde a la proyeccion en la superficie del piano de falla ILB. 
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Espectros de respuesta para la Zona 3 



[l+a (CB13)-FW 
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Figura 3 Espectros encontrados para 84% de nivel de confiabilidad para cada una 
de las micro zonas del Centra Norte aplicando los modelos de C & B (2013); 
ASK (2013), Zhao etal. (2006) y espectro ponderado. 



La microzona 1 es la de mayor peligrosidad sismica debido a que se halla 
directamente sobre la falla. La zona 4 es la que tiene las menores ordenadas 
espectrales debido a que esta en el Footwall. 

La informacion del tipo de suelo se obtuvo del estudio de microzonificacion 
sismica de Quito por ERN (2012), Aguiar (2013) Se tienen estudios de suelo en 
varios puntos de la ciudad de donde se obtuvo la velocidad de la onda de corte en 
los 30 primeros metros, V s30 y en los puntos donde no se tenfan valores se 
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interpolo linealmente con los sitios donde hay estudios de suelo. El efecto de sitio 
es fundamental considerar en la determinacion de las formas espectrales. 

Tabla 2 Espectros para las cuatro micro zonas del 
Centra Norte de Quito 



T(seg) 


Aceleracion es 


Declral (m/seg 2 ) 


Zona 1 


Zona 2 


Zona 3 


Zona 4 


PGA 


7.88 


7.40 


7.09 


5.95 


0.05 


9.96 


9.40 


9.01 


7.60 


0.1 


14.79 


14.01 


13.42 


1 1 .62 


0.15 


16.80 


15.89 


15.42 


13.32 


0.2 


17.77 


16.78 


16.10 


14.01 


0.25 


18.31 


17.16 


16.39 


14.06 


0.3 


18.60 


17.35 


16.27 


13.93 


0.35 


17.81 


16.55 


15.38 


13.12 


0.4 


16.42 


15.20 


14.13 


1 1 .96 


0.45 


15.14 


14.00 


13.03 


10.95 


0.5 


14.06 


12.99 


12.1 1 


10.1 1 


0.55 


13.07 


12.08 


1 1 .28 


9.35 


0.6 


12.17 


1 1 .23 


10.52 


8.66 


0.65 


1 1 .34 


10.46 


9.84 


8.03 


0.7 


10.58 


9.75 


9.21 


7.46 


0.75 


9.89 


9.09 


8.65 


6.93 


0.8 


9.25 


8.49 


8.13 


6.46 


0.85 


8.67 


7.94 


7.66 


6.03 


0.9 


8.14 


7.44 


7.22 


5.63 


0.95 


7.64 


6.97 


6.82 


5.27 


1 


7.19 


6.54 


6.44 


4.94 


1 .05 


6.76 


6.13 


6.08 


4.63 


1.1 


6.36 


5.76 


5.74 


4.34 


1.15 


5.99 


5.41 


5.42 


4.08 


1 .2 


5.65 


5.09 


5.12 


3.84 


1.25 


5.33 


4.80 


4.84 


3.62 


1 .3 


5.04 


4.52 


4.57 


3.41 


1.35 


4.76 


4.27 


4.32 


3.22 


1 .4 


4.50 


4.03 


4.08 


3.04 


1.45 


4.27 


3.82 


3.87 


2.88 


1 .5 


4.05 


3.61 


3.67 


2.73 


1.55 


3.84 


3.43 


3.48 


2.60 


1.6 


3.65 


3.26 


3.31 


2.47 


1.65 


3.47 


3.10 


3.15 


2.35 


1.7 


3.31 


2.95 


3.00 


2.24 


1.75 


3.15 


2.82 


2.87 


2.14 


1.8 


3.01 


2.69 


2.75 


2.05 


1.85 


2.88 


2.58 


2.64 


1.96 


1.9 


2.76 


2.47 


2.54 


1.89 


1.95 


2.65 


2.38 


2.45 


1.81 


2 


2.55 


2.29 


2.37 


1.75 


2.05 


2.46 


2.21 


2.29 


1.69 


2.1 


2.37 


2.14 


2.22 


1.64 


2.15 


2.29 


2.07 


2.16 


1.58 
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2.2 


2.22 


2.01 


2.10 


1.54 


2.25 


2.15 


1.96 


2.05 


1.50 


2.3 


2.09 


1.91 


2.00 


1.46 


2.35 


2.03 


1.86 


1.96 


1.42 


2.4 


1.98 


1.81 


1.91 


1.38 


2.45 


1.93 


1.77 


1 .87 


1.35 


2.5 


1.88 


1.72 


1 .83 


1.32 


2.55 


1.83 


1.68 


1.79 


1.28 


2.6 


1.79 


1.64 


1.75 


1.25 


2.65 


1.74 


1.59 


1.70 


1.22 


2.7 


1.69 


1.55 


1.66 


1.18 


2.75 


1.65 


1.50 


1.61 


1.14 


2.8 


1.60 


1.45 


1.56 


1.10 


2.85 


1.54 


1.39 


1.50 


1.06 


2.9 


1.49 


1.33 


1.44 


1.02 


2.95 


1.43 


1.27 


1.37 


0.97 


3 


1.37 


1.20 


1.29 


0.91 1 



En la tabla 2 se presenta las ordenadas espectrales de aceleraciones para 
las cuatro microzonas del Centra Norte de Quito. 

3 ESPECTROS NEC-1 1 Y ERN-1 2 

Los espectros presentados en el apartado anterior son para un sismo 
cortical, ahora bien Quito puede ser afectada tambien por sismos interface de bajo 
angulo de buzamiento e intraplaca. Los ultimos no producen tanto dano por la 
profundidad a la que se registran pero los sismos interface si pueden generar 
dano. 

El espectro de la Norma Ecuatoriana de la Construccion NEC-1 1, al igual 
que el espectro hallado en el estudio de microzonificacion sismica para Quito, fue 
obtenido a partir de un estudio de peligrosidad sismica en forma probabilfstica 
considerando las zonas fuentes como grandes volumenes de litosfera pero con 
estos modelos se consideran las tres fuentes sfsmicas: cortical, interface e 
intraplaca. 

Por esto las estructuras situadas en el Centra Norte de Quito deben ser 
analizadas para los espectros presentados en el apartado anterior y para los 
espectros del NEC-1 1 y ERN-1 2. Concretamente para las estructuras que tienen 
un periodo fundamental corto los espectros generados en la falla ILB seran los 
mas criticos pero en cambio para las estructuras con periodos medios y altos se 
deben considerar los espectros del NEC-1 1 y ERN-12. Las ecuaciones que 
definen el espectro del NEC-1 1 son las siguientes: 

S a (g)=zF a (i + Ctj-I)^) t<T 0 

S a (g)=i 1 zF a T 0 <T<T C (2) 



S a {g)=r ] zF a g) 



T> T c 
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Donde S a (g) es la aceleracion del suelo como una fraccion de la 
gravedad; z es el factor de zona sfsmica indicado en el mapa de peligrosidad 
sfsmica del Ecuador, para Quito z = 0.4; F a ,F d ,F s , son factores de sitio indicados 
en las tablas 3 a 5. Los perfodos T 0 ,T C , que definen las ramas del espectro se 
hallan con las ecuaciones: 



T 0 = 0.1 F s 



Fa (3) 

T c = 0.55 F s 



Fa 
F a 



En las ecuaciones que definen las ramas del espectro elastico; el valor de 
r = 1, para perfiles de suelo A, B o C y r = 1.5, para perfiles de suelo D o E. Del 
analisis de las ordenadas espectrales de peligro uniforme en roca para 475 anos y 
al normalizarlos para la aceleracion maxima del suelo z, se encontro r\ = S a /z, que 
tienen los siguientes valores: r\ = 1.8 Para las Provincias de la Costa; r\ = 2.48 
Para las Provincias de la Sierra y r\ = 2.6 Para las Provincias del Oriente. 



Tabla 3 Tipo de suelo y Factores de sitio Fa 





Zona 
sfsmica 


1 


II 


III 


IV 


V 


VI 


Tipo de 
perfil del 
suelo 


Valor Z 
(Aceleracion 
esperada en 

roca, 'g) 


0.15 


0.25 


0.3 


0.35 


0.40 


>0.5 


A 


0.9 


0.9 


0.9 


0.9 


0.9 


0.9 


B 


1 


1 


1 


1 


1 


1 


C 


1.4 


1.3 


1.25 


1.23 


1.2 


1.18 


D 


1.6 


1.4 


1.3 


1.25 


1.2 


1.12 


E 


1.8 


1.4 


1.25 


1.1 


1.0 


0.85 


F 


Ver Nota 


Ver Nota 


Ver Nota 


Ver nota 


Ver Nota 


Ver Nota 



Los factores de sitio del NEC-11, son muy generales, por ese motivo es 
recomendable emplear los factores de sitio hallados para la ciudad de Quito que 
fueron hallados por ERN-12 ya que estos responden a los estudios de suelo que 
se realizaron en diferentes sitios de la ciudad, a las curvas de rigidez y 
amortiguamiento de las cangahuas de la ciudad y a estudios dinamicos no lineales 
con los que se hallaron los factores de amplificacion de los suelos de Quito, los 
mismos que fueron normalizados de tal manera que se puedan utilizar las 
ecuaciones que definen las ramas del espectro elastico del NEC-11, con la unica 
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salvedad es que ERN-12 lo hizo considerando un valor de r = 1, para cualquier 
tipo de suelo. 



Tabla 4 Tipo de suelo y Factores de sitio Fd 



Tipo de 
perfil del 
suelo 


Zona 
sismica 


1 


II 


III 


IV 


V 


VI 


Valor Z 
(Aceleracion 
esperada en 

roca, 'g) 


0.15 


0.25 


0.3 


0.35 


0.4 


>0.5 


A 


0.9 


0.9 


0.9 


0.9 


0.9 


0.9 


B 


1 


1 


1 


1 


1 


1 


C 


1.36 


1.5 


1.4 


1.35 


1.3 


1.25 


D 


1.62 


1.7 


1.6 


1.5 


1.4 


1.3 


E 


2.1 


1.75 


1.7 


1.65 


1.6 


1.5 


F 


Ver Nota 


Ver Nota 


Ver Nota 


Ver nota 


Ver Nota 


Ver nota 


Tabla 5 Tipo de suelo y Factores del comportamiento inelastico del subsuelo Fs 


Tipo de 
perfil del 
suelo 


Zona 
sismica 


I 


II 


III 


IV 


V 


VI 


Valor Z 
(Aceleracion 
esperada en 

roca, 'g) 


0.15 


0.25 


0.3 


0.35 


0.4 


>0.5 


A 


0.75 


0.75 


0.75 


0.75 


0.75 


0.75 


B 


0.75 


0.75 


0.75 


0.75 


0.75 


0.75 


C 


1.00 


1.1 


1.2 


1.25 


1.3 


1.45 


D 


1.2 


1.25 


1.3 


1.4 


1.5 


1.65 


E 


1.5 


1.6 


1.7 


1.8 


1.9 


2.0 


F 


Ver Nota 


Ver Nota 


Ver Nota 




Ver nota 


Ver Nota 



Nota.- Son suelos muy malos que requieren estudios especiales. 



En la figura 4 se presentan los factores de sitio F a ,F d ,F s hallados por ERN- 
12, el rango de variacion de estos factores difiere con relacion a los factores de 
sitio del NEC-11. Evidentemente que es una mejor opcion utilizar los factores de 
sitio de ERN-12 ya que como se indico anteriormente responden a las condiciones 
locales de los suelos de Quito. 
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Figura 4 Factores de sitio F a ,F d ,F s hallados para la ciudad de Quito por ERN-12 

Ahora bien en las figuras 5 y 6 se presentan los espectros hallados para 
las zonas 1 y 2, para un sismo de magnitud 6.2 en el segmento de falla ILB 
(Espectro Ponderado) y se presentan ademas los espectros que se hallan en 
diferentes sectores del Centra Norte de Quito, utilizando los factores de sitio 
propuestos por ERN-12. 

En las figuras 5 y 6 se aprecia que las ordenadas del espectro ponderado 
son mayores, aproximadamente en un 30% con respecto a las ordenadas 
espectrales que se hallan con ERN-12, pero esto es unicamente en la rama de 
penodos cortos, menores a 0.5 seg. 

Para la zona de penodos intermedios y grandes la situacion es diferente 
las ordenadas espectrales que se hallan con ERN-12 son mas altas. Esto es muy 
importante tener en cuenta para que no se piense que lo mas seguro es utilizar el 
espectro de ILB para el diseho de cualquier tipo de estructuras del centra norte de 
Quito por que se tienen ordenadas espectrales mucho mas grandes en la zona de 
penodos cortos. 

La Normativa Sismica de Chile de 1996 responde a sismos corticales y los 
edificios que sufrieron gran daho en el mega sismo del 2010 fueron disehados con 
esta normativa y tuvieron gran daho porque los espectros del sismo en la zona de 
penodos intermedios y largos tenian varios picos caracteristicos cuyas 
aceleraciones eran notablemente mayores a las de la Norma de 1996. (Moroni et 
al. 2012, Aguiar 2010). 
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Espectros de respuesta para laZona 1 




Periodo(s) 

Figura 5 Espectro ponderado encontrado en el estudio para la Zona 1 y los 
Espectros que se hallan al aplicar los factores de sitio de la 
Microzonificacion sfsmica de Quito. (ERN-12) 

Espectros de respuesta para la Zona 2 











3 onderado 
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ERN12- Voz de los Andes 
475anos) 










ERN12- Chaupicruz (475anos) 
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0 0.5 1 1.5 2 2.5 3 



Periodo (s) 

Figura 6 Espectro ponderado encontrado en el estudio para la Zona 2 y los 
Espectros que se hallan al aplicar los factores de sitio de la 
Microzonificacion sfsmica de Quito. (ERN-12) 



Con este artfculo lo que se pretende es que los Proyectistas estructurales 
consideren el espectro de ILB para el diseno de edificios de poca altura, cuyo 
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periodo fundamental sea menor a 0.5 seg., y para edificios de mayor altura lo 
mejor es utilizar el espectro de ERN-12. 

4 ANALISIS ESTATICO NO LINEAL 

El analisis estatico no lineal, tambien conocido como la Tecnica del 
Pushover o la Tecnica del Empujon, aparecio en la decada de los anos 90 del 
siglo pasado y en el presente siglo ha habido importantes contribuciones 
tendientes a su optimizacion entre las que se destacan las propuestas por Chopra 
y Goel (2001, 2002), Goel y Chopra (2004). Sin embargo de ello un buen numero 
de profesionales solo saben utilizar programas como el DRAIN-2DX para 
encontrar la curva de capacidad sismica resistente que relaciona el cortante basal 
V, con el desplazamiento lateral maximo D t empleando analisis estatico no lineal; 
por este motivo en este apartado se presentan aspectos generales de la forma 
como se halla esta curva, orientados al uso del sistema de computacion CEINCI- 
LAB. 

4.1 Modelo Bilineal del diagrama M — 0 

La no linealidad de los elementos viene dado por el diagrama momento- 
curvatura, (M - 0); en este artfculo se trabaja con un modelo bilineal definido por 
los puntos Y, U; y esta orientado a estructuras de acero. Donde el punto Y es el 
punto de fluencia y U el punto de resistencia ultima o de plastificacion. De tal 
manera que solo se consideran dos rigideces a flexion, una para el rango elastico 
(EI) e y otra para el rango plastico (EI) p , verfigura 7. 

(EI)e = / 

0Y (4) 

P 0y -0 Y 

Donde M Y ,0 Y son el momenta y curvatura a nivel de fluencia; M Ul 0 u son el 
momenta y curvatura a nivel de plastificacion. En la figura 7 se aprecian dos 
diagramas uno que se halla en el primer cuadrante, para cuando la parte inferior 
de la seccion transversal trabaja a flexion y otro en el tercer cuadrante para 
cuando la parte superior de la seccion transversal trabaja a flexion. En elementos 
con seccion transversal simetrica los dos diagramas (M - 0) son iguales; en este 
artfculo se trabaja el caso de secciones simetricas. 

Para cada uno de los elementos se debe indicar el momenta y curvatura 
de los puntos Y, U, en el nudo inicial (a), centra de luz (o) y nudo final (b), ademas 
de eso se da como data la rigidez axial EA, de tal manera que se dan trace datos 
para cada uno de los elementos en un archivo denominado Tablal. En la tabla 6 
se presenta el contenido de este archivo de datos para ser utilizado en el sistema 
de computacion CEINCI-LAB. 
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Figura 7 Modelo bilineal del diagrama momento curvatura. 



Tabla 6 Datos del diagrama momento curvatura para un elemento. 



Nudo Inicial (a) 


Centra de Luz (o) 


Nudo final (b) 


Axial 


My 


My 


0 Y 


0(7 


M Y 




0 Y 


0u 


M Y 




0 Y 


0u 


EA 



4.2 Modelo de plasticidad concentrada 

Se trabajo con el modelo de plasticidad concentrada de Giberson (1969) 
que concentra la rigidez a flexion en tres puntos como se aprecia en la figura 8. 
Siendo (£7) a para el nudo inicial; (£7) 0 para el centra de luz y (£7) b para el nudo 
final. La matriz de flexibilidad para el modelo indicado en la figura 8, en que no se 
consideran los desplazamientos como cuerpo ngido es la siguiente. 

(1 + S a ) -1 
-1 (1 + S b ). 

(5) 

, (W)o (£Oo 

a ma b m\ 



u(EI) a (El) 0 (El) b 




L 



Figura 8 Modelo de Plasticidad Concentrada de Giberson 

Posteriormente se incorpora el efecto axial al modelo de la figura 8, con lo 
que la matriz de flexibilidad es de 3 por 3, la inversa de esta matriz es la matriz de 
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rigidez k v En Aguiar (2014) a ese sistema se conoce como sistema 1; el sistema 2 
es en coordenadas locales, cuya matriz de rigidez k 2 = T\_ 2 * fe x * T ± _ 2 y en el 
sistema 3 en coordenadas globales, la matriz de rigidez fe 3 = r|_ 3 * ^2 * T 2 _ 3 . Las 
matrices T^_ 2 y T 2 -3 son l as matrices de paso de los sistemas 1 al 2 y del 2 al 3. A 
esta ultima tambien se le conoce con el nombre de matriz de rotacion. El programa 
que determina la matriz de rigidez en el sistema 3 a partir de la matriz de 
flexibilidad en sistema 1 se denomina kmiembro_giberson_po 

Cuando el nudo inicial o final se encuentran en el rango elastico la rigidez 
a flexion es (£7) e y es valida cuando el momenta actuante es menor o igual al 
momenta de fluencia. Cuando se supera este punto la rigidez a flexion vale (£7) p , 
y es valida hasta cuando el momenta actuante es menor a M u . Para cualquiera de 
estos dos casos se trabaja con el programa kmiembro_giberson_po. 

Ahora bien, cuando la seccion Mega al momenta M 0 se forma una rotula 
plastica, en este caso la seccion ya no es capaz de absorber mas momenta y solo 
puede rotar. 



(El)- 



(El) b 



II (ED, 



(El) c 



WWW fife 



Rotula Plastica en Nudo Inicial 



Rotula Plastica en Nudo Final 



Figura 9 Modelos para cuando se forma rotula plastica en el nudo inicial o final, en 

el sistema 1 

A la izquierda de la figura 9, se presenta el caso en que la rotula plastica 
se forma en el nudo inicial, en este caso solo se admite incremento de momenta a 
flexion en el nudo final. La matriz de flexibilidad para el sistema 1, es: 



/ = 



6 (£/)„ 



[1 + s b ] 



(6) 



El programa que determina la matriz de rigidez de un elemento en 
coordenadas globales, cuando se forma una rotula plastica en el nudo inicial, es: 
kmiembro_giberson_po_sinSa. 

Cuando la rotula plastica se forma en el nudo final, el modelo de analisis 
es el indicado a la derecha de la figura 9. La matriz de flexibilidad en sistema 1 , es: 



6 (£/)„ 



[1 + S«] 



(7) 



En este caso el programa es: kmiembro_giberson_po_sinSb. Se pudo 
trabajar con una matriz de dos por dos para cuando se forma rotula plastica en el 



Roberto Aguiar, David Mora y Enrique Morales 



18 



nudo inicial o en el nudo final. Para el primer caso la primera fila y columna tendra 
solo ceros y el termino (2,2) la cantidad indicada en la ecuacion (6). Para el caso 
de que la rotula plastica se forme en el nudo final, la cantidad indicada en la 
ecuacion (7) va en la posicion (1,1) y la segunda fila y columna tiene ceros. 

4.3 Incrementos de carga 

La curva de capacidad sfsmica resistente depende de la forma como se 
aplica la carga lateral en cada uno de los pisos. Lo mas adecuado es que esta 
carga sea proporcional al primer modo de vibracion 0 A ; debido a que la mayor 
parte de estructuras regulares trabajan en el primer modo de vibracion. 

Otra forma de hacerlo, es encontrar las cargas en cada piso proporcional 
al primer modo y al peso de cada piso w t . De tal manera que la fuerza lateral en el 
piso Fj se halla con la siguiente ecuacion, propuesta por el ATC-40 

Donde V 0 es el cortante basal que uno se impone para cada uno de los 
ciclos de carga. Se destaca que se esta trabajando solo con el primer modo de 
vibracion, de tal manera que el subfndice t de la ecuacion (8) hace referencia al 
piso. Se puede encontrar la curva de capacidad sfsmica resistente asociada al 
segundo modo o al tercer modo, trabajando con el segundo y tercer modo. Chopra 
yGoel(2001). 

4.4 Diagrama de flujo de curva de capacidad sfsmica resistente 

En la figura 10 se presenta el diagrama de flujo de la curva de capacidad 
sfsmica resistente; en la primera parte se dan los datos para la generacion del 
portico utilizando los programas de CEINCI-LAB para porticos pianos. Se inicia 
obteniendo la matriz que contiene a los grados de libertad de la estructura, 
denominada CG; posteriormente se determinan dos vectores que contienen los 
nudos iniciales y finales del portico, Nl, NJ; se da las coordenadas de los nudos y 
se encuentran los vectores X, Y. Luego se obtienen las propiedades geometricas 
de los elementos y los vectores de colocacion; esta es la informacion basica para 
empezar el analisis. Se destaca que no se considera el efecto P — A. 

Posteriormente en un gran lazo que empieza con _/' > 0 se describe el 
analisis estatico incremental, el mismo que finaliza cuando la deriva global de la 
estructura es mayor a la cantidad prefijada, que es funcion de la altura total del 
edificio o cuando la matriz de rigidez es singular. 
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PROGRAMA Pushover C INIpIO J 

Tablal-Matriz de momento curvatura en inicio, medio y final de c/elem 
nod-Num de nudos; np-Num de pisos; nr-Num de nudos empotrados 
GEN-[i,ia,ib,nig,ii,ina,inb];NUDOS-[i,xi,yi,ij,inci,dx,dy];Fm-[Num de miembro cargado,Carga repartida] 



progrCEINCI-LAB:[CG,ngl]=cg_sismo(nod,np,nr); 
CG -Matriz de coordenadas generalizadas 
ngl-Numero de grados de libertad 

prgr CEINCI-LAB:[L,seno,coseno]=longitud (X,Y,NI,NJ); 

L=Vector con longitudes de los elementos 
seno-Vector con sen del angulo de inclinacion de c/elem 
coseno-Vector con cos del angulo de inclinacion de c/elem 



prgrCEINCI-LAB:[NI,NJ]=gn_portico(GEN); 
GEN-[i ia ib nig ii ina inb] 
NI,NJ-Vectores de nudos inic. y finales 



prgrCEINCI-LAB:[X,Y]=glinea_portico(NUDOS) 
X,Y-Vector con coordenadas de los nudos 



prgrCEINCI-LAB:[VC]=vc(NI,NJ,CG); 
VC-Matriz con vectores horiz de colocacion por elem 



mbr-tamano de vector L;mbr-Numero de elementos 
DeltaQ-0.5; Increm. de Cort. Basal en c/iteracion 



Ht- maximo valor de vetor Y; 

prgr CEINCI-LAB:[Mt,M]-masaa(Fm,NI,Y,L,np} 
Mt-lvlasa total de estructura ;M=Matriz de masas de estructura 



Tabla2- Matriz de ceros de mbr x 6; Tabla3-Matriz de ceros de mbr x 8 
icon=iconu-derg-conty-cont1My-j-VV-Dtt=0 
Qa-Vector vert de ceros de np filas 



HZ 

j>=0 

Comienza bucle, se detiene cuando la 
estructura colapsa 



b>| j=i + i 



progrCEINCI-LAB: 

[K,icon,iconu,Tabla3,VC,ngl]-krigidez_giberson_pofngl,Tabla1 ,Tabla2,Tabla3,L,seno,coseno,VC,E,icon,iconuj) 
K-Matriz de rigidez de la estructura;Tabla2=Contine N,V,M,N , ,V , ,M' de cada elemento 

Tabla3-Col. 1 instante en que pasa a My en el nudo inicial; 
Col 2 inst en que pasa a My el nudo final;Col 3 inst en que pasa a Mu el nudo inicial 
Col 4 inst en que pasa a Mu el nudo final;Columna 5 secuenc de dano a My en nudo inicial 
Col 6 secuenc de dano a My en nudo final;Col 7 secuenc de dano a Mu en nudo inicial. 
Col 8 secuenc de dano a Mu en nudo final; icon-Contador de fluencia de un elemento. 
iconu-Contador de formacion de rotula plastica de un elem;j-ldentifica ciclo de carga. 




progrCEINC 

[Q,V,Q2]-cargas_po(Qa,ngl,Fm,L,seno,coseno,VC); 

Qa=Vector vert de cargas horizontales en c/piso 
Q-Vector de cargas generalizadas;V-Cortante Basal 
Q2-Matriz con acciones de empotramiento perfecto 



Desplaz. y giros de grados de libertad 
q=K\Q; 



Desplazamiento en el tope 
Dt-valor de vector q en GDL np; 



progr UbINCI-LAB: 

[Sdd,Saa,fi1]=espe_cap(K,M,Mt,np,ngl,Dtt,W); 
Sd(j)=Sdd;SaG)=Saa;EspectroC(j,1)=Sd(j);EspectroC(j,2)=Sa(j) 



derg=Dt/Ht(Deriva global);VV^VV+V(V Basal acumulada) 
Dtt-Dtt+Dt;Desplaz en tope de cada iteracion 



Colapso por deriva 
global>=0.02 



progr CEINCI-l_AB:[FF(:,:j)]^fuerzas_po(ngl,Tabla1,Tabla2,Tabla3,L,seno,coseno,VC,q,Q2) 
F-Matriz tridimensional de Fuerzas y Momentos finales en cada elemento en coordenadas locales acumulada en } 



I Tabla2=(FF de iteracion j)+Tabla2 





contl My-cont1 My+1 ; 
DeltaQ=DeltaQ/100 
Tabla2-Tabla2-FF(:,:j-1) 
Tabla3-Matriz de ceros mbrx8 
iconu-Q;i-j-2 



contl My=3 — e> 



contl My-cont1 My+1 ; 
DeltaQ=DeltaQ*100 
V y Despl de Pto. de My: 

Vacy=Vac(j)=Vac(i-1 ) 
Dtacy^Dtac{j)^Dtac(j-1 ) 



Qa(:,j+1)=(sum(M)'.*fi1*DeltaQ)/sum(sum(M)'.*1i1) 



-I— 



Figura 10 Diagrama de flujo de la curva de capacidad sismica resistente 
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En el primer lazo donde se llama al programa krigidez_giberson_po ahf 

se indica la forma como se van numerando, en primer lugar las secciones que 
superan el momento de fluencia y posteriormente el momenta plastico, informacion 
que es almacenada en Tabla3. Se destaca que este lazo se encuentra dentro de 
un gran lazo. Por lo tanto corresponde a un determinado ciclo de carga; conforme 
se incrementa la carga otras secciones ingresaran al rango no lineal. 

Posteriormente se indica la forma de calculo del vector de cargas 
generalizadas y la solucion del portico mediante el Metodo de los 
Desplazamientos, destacando que en el primer ciclo de carga se considera 
unicamente la carga vertical que gravita en los porticos. A partir del segundo ciclo 
de carga se aplican las cargas laterales. 

En la parte final del lazo a la derecha se presenta la forma en que se 
reduce el incrementa de carga, cuando se pasa de la rama elastica a la inelastica; 
se reduce el incrementa de carga para llegar en forma exacta al punto M Y . Luego 
de esto se vuelve al incrementa de carga con que se venfa calculando. 

En la figura 11 se presenta el calculo de la matriz de rigidez por 
ensamblaje directo, para ello en base a los momentos actuantes se pregunta si 
superaron el momento de fluencia, en caso afirmativo se cambia la rigidez a 
flexion del rango elastico al rango plastico, en el nudo del elemento que se esta 
analizando. Esto se lo hace en cada uno de los elementos tanto en el nudo inicial 
como final. 

Cuando se cambia de rama en el modelo bilineal, se regresa al punto 
antes del cambio de rama y se disminuye el incrementa de carga lateral para llegar 
en forma exacta al punto de cambio de rama. De tal manera que en cada ciclo de 
carga se actualiza la rigidez a flexion del nudo inicial, centra de luz y nudo final; 
procediendo a calcular la matriz de rigidez del elemento en coordenadas globales 
y ensamblar la matriz de rigidez de la estructura. 

Cuando una seccion de un elemento se encuentra en el rango plastico la 
pregunta que se hace es: el momento actuante en dicho punto supero al momento 
M u en este caso se formo la rotula plastica. Esta parte del programa se encuentra 
en la mitad del diagrama de flujo indicado en la figura 11 y se aprecia que cuando 
se forma la rotula plastica en el nudo inicial se llama al programa: 
kmiembro_giberson_po_sinSa. Por otra parte, cuando la rotula plastica se forma 
en el nudofina se llama al programa: kmiembro_giberson_po_sinSb. Finalmente 
cuando se tiene un elemento con dos rotulas plasticas, en el nudo inicial y final, se 
llama al programa: kmiembro_giberson_po_sinSaSb. 

En la parte final de la figura 11 se presenta el ensamblaje directo de la 
matriz de rigidez y para ello se trabaja con el vector de colocacion VC de cada 
elemento. En el diagrama de flujo se observa que se ensambla toda la matriz de 
rigidez y se trabaja con un arreglo de dos dimensiones. Una optimizacion seria 
ensamblar solo la matriz triangular superior y trabajar con un arreglo lineal la 
matriz de rigidez. En Aguiar (2014) se presentan varias opciones de ensamblaje 
de la matriz de rigidez, considerando ancho de banda constante o ancho de banda 
variable (Tecnica del Skyline) esto es adecuado para cuando se tienen estructuras 
muy grandes. 
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PROGRAMA krigidez_giberson_po ( INICIO 



Tablal-Matriz de momento curvatura en inicio, medio y final de cada elemento por fila;Tabla2-Contine N,V,M,N',V',M' de cada elemento por fila 
Tabla3-c/fila representa un elem (Columnas descritas en progr PUSHOVER;ngl-Numero de grados de libertad; L-Vector con longitudes de los elem . 
seno-Vector con seno del angulo de inclinacibn de cada elemento;coseno-Vector con coseno del angulo de inclinacion de cada elemento 
icon-Contador de fluencia de algiin elemento. ;iconu-Contador de formacibn de rotula plastica de algun elemento. 
VC-Matriz con vectores horizontal de colocacion por elementoj-ldentifica el ciclo de carga 



mbr-tamano de vector L; 
mbr-Numero de elementos 



SS-Matriz cuadrada de ceros de tamano de los gradoslde 
libertad 



>- i=1,mbr ^ >> 

, if 

IMya-Tabla3fi,5);iMyb-Tabla3(i,6);iMua-Tabla3fi,7) 
iMub=Tabla3(i,8);M=[Tabla2(i,3)|->M en nudo inicial 
Mp=|Tabla2(i,6)|-->M en nudo final 



M<=Tabla1(i,1) 

Tabla3(i,1 )=0^> 

& Tabla3( ,3)-0 

Jsi 

Ela-Tabla1(i,1)/Tabla1(i,3) 



M>Tabla1(i,2)^ N0 
3 Tabla3{i,3)<- " 



E , (Tabla1(i,2i-Tabla1(i,1j) 
(Tabla1(i.3VTabla1(i.3» 





_ ;j:SI 

icon=icon+1 
Tabla3(i,5)=icon;Tabla3(i,1Hj 




iconu-iconu+1 
Tabla3(i,7)=iconu;Tabla3(i,3)=jj 




Elo-Tabla1(i,5)/Tabla1(i,7) 
EA=Tabla1(1,13) 
EIF-[Ela Elo Elb]-->Vector con rigideces de elem 
lon-l_(i);sen-seno(i);cosen-(i); 




Figura 11 Determinacion de la matriz de rigidez de la estructura por 
ensamblaje directo. 
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PROGRAMA kmiembro giberson INICIOJ 



- EIF=[Ela Elo Elb]-->Vector con rigideces a flexion de elem 

EA->Rigideces axial de elem 
L=Longitud del elemento;seno=Sen del angulo de inclinacion .. 
del elem;coseno=Cos del angulo de inclinacion del elem . 

I 



Ela=|EIF<1,1)|; Elo=|EIF(1,2)|; Elb=|EIF{1 ,3)| 



Sa=Elo/Ela 
Sb=Elo/Elb 



flex=L/(6*Eolo]*[(1+Sa) -1 
-1 (1+Sb) 



T12=[0 1/L 1 0 -1/L 0; 
0 1/L 0 0 -1/L 1 
-1 0 0 1 0 0] 



K1=T12?*K1*T12 



T23= [ 








0 0 0 0 










0 0 0 0 




0 


0 




10 0 0 




0 


0 


0 


coseno seno 0 




0 


0 


0 


-seno coseno 0 




0 


0 


0 


0 0 1] 



PROGRAMA kmiembro_giberson_po_sinSaSb INICIO 



/ EIF=[Ela Elo Elb]~>Vector con rigideces a flexion de elem 

EA~> Rigideces axial de elem 
L=Longitud del elemento;seno=Sen del angulo de inclinacion , 
del elem;coseno=Cos del angulo de inclinacion del elem/ 



T12=[-1 0 0 1 0 0] 



123= [coseno 


seno 


0 0 


0 


0 


-seno cc 




0 0 


0 


0 


0 


0 


1 0 


0 


0 


0 


0 0 






0 


0 


0 0 






0 


0 


0 0 


0 0 


1] 





K3=T23>K2'T23 



PROGRAMA 
kmiembrogibersonposinSa 



PROGRAMA 
kmiembrogibersonposinSb 



INI CIO 



EIF=[Ela Elo Elb]->Vector con rigideces a flexion de elem 
EA«>Rigideces axial de elem 
L=Longitud del elemento;seno=Sen del angulo de inclinacibn 
del elem;coseno=Cos del angulo de inclinacibn del elem 



Elo=|EIF(1,2)| 
Elb=|EIF(1,3)| 

i 

Sb=Elo/Elb 

3E 



flex=L/(6*Eolo)*[(1+Sb) 




EIF=[Ela Elo Elb]->Vector con rigideces a flexion de elem 
EA->Rigideces axial de elem 
L=Longitud del elemento;seno=Sen del angulo de inclinacion , 
del elem;coseno=Cos del angulo de inclinacion del elem/ 

EIO=|EIF(1,2)| 
Ela=|EIF(1,3)| 

I 



flex=L/(6*Eolo)*[(1+Sa) 



T12=[0 1/L 0 0 -1/L1 
-1 0 0 1 0 0] 



T23= [coseno seno 
-seno coseno 



10 0 0 
0 coseno seno 0 
0 -seno coseno 0 
0 0 0 1] 



T12=[0 1/L 1 0 -1/L0 
-1 0 0 1 0 0] 



K1=T12>K1*T12 



T23- [coseno 




o 0 0 0 0 






o 0 0 0 0 


0 


0 


10 0 0 


0 


0 


0 coseno seno 0 


0 


0 


0 -seno coseno 0 


0 


0 


0 0 0 1] 



Figura 12 Determinacion de la matriz de rigidez de un elemento. 
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En la figura 12 se presentan cuatro diagramas de flujo para el calculo de la 
matriz de rigidez de un elemento. El primer diagrama es para cuando en el 
elemento no se han formado rotulas plasticas; el segundo cuando se han formado 
rotulas en el nudo inicial y final; el tercero, a la izquierda cuando la rotula se forma 
en el nudo inicial y el cuarto, a la derecha cuando la rotula se forma a la derecha. 



5 ESTRUCTURA DE ANALISIS Y CAPACIDAD SISMICA 



Se desea encontrar la curva de capacidad sismica resistente que relaciona 
el cortante basal con el desplazamiento lateral maximo de la estructura de acero 
de 6 pisos indicada en la figura 13. Si el criterio de colapso esta definido por la 
deriva global y g expresada como un porcentaje de la altura total del edificio H T . Se 
halla esta curva para: i) y g = 0.02ff r ;ii) y g = 0.03 H T ; iii) y g = 0.04 H T ; iv) y g = 
0.05 Hr. 



r 



• — i 

x. 



i 



i 

X 
5 



rH 

o 

i 



W24 X 76 



W24 X 76 



W27 X 94 



-4- W27X94 



W30 X 99 



W30 X 99 



3 @ 7.315 m = 21.945 m 



£ 



1.852 T/m 
3.269 T/m 

3.269 T/m 
3.269 T/m 

3.269 T/m 
3.269 T/m 



Figura 13 Estructura de seis pisos y carga uniforme repartida en cada piso. 
Tabla 7 Diagrama Momento Curvatura y rigidez axial para las columnas 



Perfil 


P 


M Y 


My 


0y 


0u 


EA 




(T.) 


(Tm.) 


(Tm.) 


(1/m) 


(1/m) 


(T.) 


W14X109 


599 


91.20 


109.44 


0.00884 


0.0972 


432600 


W24X104 


573 


137.00 


164.40 


0.00532 


0.0585 


413700 


W14X159 


874 


136.00 


163.20 


0.00862 


0.0948 


632100 


W27X146 


803 


219.00 


262.80 


0.00468 


0.0514 


581700 


W14X193 


1060 


169.00 


202.80 


0.00844 


0.0929 


768600 


W30X173 


951 


288.00 


345.60 


0.00421 


0.0463 


688800 
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En la tabla 7 se presenta el diagrama momento curvatura para los perfiles 
tipo de las columnas, para la carga axial indicada en la segunda columna; en la 
ultima columna se indica la rigidez axial. En la tabla 8 se tiene lo mismo pero para 
las vigas, en este caso el diagrama M - 0 no depende de la carga axial. 

Tabla 8 Diagrama Momento Curvatura y rigidez axial para vigas 



Perfil 


My 


M u 


0y 


0u 


EA 




(Tm.) 


(Tm.) 


(1/m) 


(1/m) 


(T.) 


W24X76 


95.00 


114.00 


0.00543 


0.0598 


304500 


W27X94 


132.00 


158.40 


0.00485 


0.0534 


375900 


W30X99 


148.00 


177.60 


0.00446 


0.0491 


394800 



2% Curva V-Deformacion terraza 




0.2 0.3 
Dt(m.) 



Esquema de Estructura y de Rotulas 
25 



20 



15 



10 



Tfl ~i~ < 



10 15 20 



Figura 14 Diagrama Momento Curvatura para criterio de colapso y g = 0.02 H T 

En la figura 14 se presenta, a la izquierda, la curva de capacidad sismica 
resistente V — D t cuando el criterio de colapso es y g = 0.02 H T , A la derecha se 
aprecia que no se formo ninguna rotula plastica. El sistema de computacion 
CEINCI-LAB representa con un cuadrado las rotulas plasticas. Varias secciones 
del portico superaron el momento de fluencia por este motivo es la curva que 
existe en V - D t . 

Cuando el criterio de colapso y g = 0.03 H T , se forman dos rotulas plasticas 
en los nudos iniciales de las columnas inferiores, ver figura 15. 



Se habfa indicado que la Tablal contiene los momentos y curvaturas de 
los puntos Y, U, en el nudo inicial, centra de luz y nudo final; ademas de ello se 
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dar como dato para cada elemento la rigidez axial. Esta tabla es suministrada por 
el usuario del sistema CEINCI-LAB. Ahora bien el programa crea y reporta la 
Tabla2 que contiene la fuerza axial, fuerza de corte y momento del nudo inicial y 
final, de cada uno de los elementos. 



3 % Curva V-Deformacidn terraza Esquema de Estructura y de Rdtulas 




Dt(m) 

Figura 15 Diagrama Momento Curvatura para criterio de colapso y g = 0.03 H T 



Tambien se crea la Tabla3 que basicamente va indicando la secuencia de 
las secciones que superan el momento de fluencia. Para primera seccion que 
supera el momento M y se determina el desplazamiento en el tope y este punto 
representa el fin del rango elastico de la estructura e inicio del rango plastico. En 
Tabla3 tambien se tiene la secuencia de la formacion de rotulas plasticas. Ademas 
se indica para los dos casos (superacion de My.My) el cortante basal. 

En la figura 16 se presenta la curva V -D t , cuando el criterio de colapso 
viene dado por y g = 0.04 H T . A la derecha de esta figura se observa que a mas de 
las rotulas plasticas de los nudos iniciales de las columnas centrales, se han 
formado rotulas plasticas en las tres vigas del primer piso, en el nudo inicial y final. 

Como se podra deducir cuando se desea obtener la curva V - D t , se debe 
definir el criterio de colapso, una buena parte de Proyectistas Estructurales 
consideran y g = 0.02 H T . Ademas de ello, en cada ciclo de incremento de carga se 
encuentra el determinante de la matriz de rigidez y si es menor a cero, la 
estructura es inestable y colapsa, de tal manera que es otro criterio de colapso que 
considera CEINCI-LAB. 
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4% Curva V-Deformacion terraza 

rT=r 




Esquema de Estructura y de Rotulas 
25 
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Figura 16 Diagrama Momento Curvatura para criterio de colapso y g = 0.04 H T 
5% Curva V-Deformacion terraza 




Esquema de Estructura y de Rotulas 
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Figura 17 Diagrama Momento Curvatura para criterio de colapso y g = 0.05 H T 



27 



Peligrosidad sismica de Quito y el Metodo del Espectro de Capacidad 



Finalmente hay otro criterio de colapso que el Proyectista debe darse 
cuenta si se ha formado un mecanismo y por lo tanto la estructura es inestable. 
Hasta ahora no se tiene ningun mecanismo de colapso; ni en la figura 17 en que el 
criterio de colapso viene dado por y g = 0.05 H T . 

En la figura 18 se presenta a la izquierda las secciones donde se formaron 
rotulas plasticas y se aprecia que todos los nudos iniciales de las columnas se han 
rotulado en consecuencia se genera un mecanismo y la estructura colapso. A la 
derecha de la figura 18, en la curva V — D t se han sehalado los instantes en que 
se formaron las diferentes rotulas plasticas. En este caso la deriva global es mayor 
al 5% de la altura global. Por cierto, la deriva global relaciona el desplazamiento 
lateral en el tope del edificio con la altura total del mismo. 



Cufva V-Deformacion terraza 
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Figura 18 Mecanismo de colapso de portico de seis pisos de acero. 

En la figura 19 se indican las 10 primeras secciones que superaron el 
momento de fluencia M Y . Notese que la estructura ingreso al rango no lineal 
cuando el cortante basal es de 212 T. 
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Curva V-Deformacion terraza 
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Figura 19 Diez primeras secciones que ingresan al rango no lineal, superan el 

Momenta de Fluencia 

6 ESPECTRO DE CAPACIDAD 



El concepto de espectro esta asociado a un sistema de un grado de 
libertad. Esto es muy importante tener en cuenta para cuando se aplique el 
Metodo del Espectro de Capacidad y se encuentre el Punto de Demanda. 

En Aguiar (2002, 2003) esta bien detallado el marco teorico para pasar de 
la curva de capacidad sismica resistente V - D t al espectro de capacidad S a - S d . 
Las ecuaciones de calculo son las siguientes. 



S . = 



Vi 



a-, M T 



Sdi — 



FPi 



(9) 
(10) 



Donde D ti ,Vi son las coordenadas de un punto de la curva de capacidad 
sismica resistente para los cuales se determina el desplazamiento y aceleracion 
espectral S di ,S ai ; M T es la masa total de la estructura; a t es el coeficiente de masa 
modal del primer modo; FP 1 es el factor de participacion del primer modo. 



«i = 



(11) 



29 



Peligrosidad sismica de Quito y el Metodo del Espectro de Capacidad 



FP.=^- (12) 
1 01 M 0i 

Las variables todavfa no indicadas son: <p 1 es el primer modo de vibracion 
normalizado a la unidad en el tope; M es la matriz de masas; J es el vector de 
incidencia del movimiento del suelo con las coordenadas laterales del portico; para 
el caso piano J es un vector unitario. Se destaca que el factor de participacion 
modal se obtiene en valor absoluto. 

Para la estructura de acero de seis pisos, mostrada en la figura 13, se 
encontro el espectro de capacidad, considerando la curva de capacidad sismica 
resistente indicada en la figura 15, que corresponde a: y g = 0.03 H T . 

En la parte central de la figura 20 se presenta el espectro de capacidad, la 
misma que fue obtenida con el programa espe_cap. Ingresa como dato, la matriz 
de rigidez de la estructura K para el punto de interes D ti ,Vi, el programa condensa 
la matriz de rigidez a las coordenadas laterales y con la matriz de masas M se 
halla los penodos y modos de vibracion; interesa solo el valor de <$> x ; los factores 
a lf FPj ; y los respectivos puntos del espectro de capacidad S ai - S di . 



3% Curva V-Deformacion tope 




Espectro de Capacidad Estructura y de Rdtulas 

25 




20 



15 



10 



_m 



10 



20 



Dt (m.) — Sd (m.) 

Figura 20 Curvas V - D t ; espectro de capacidad y formacion de rotulas plasticas 
para criterio de colapso y g = 0.03 H T . Estructura de acero de 6 pisos. 
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Para encontrar las fuerzas y desplazamientos laterales en cada piso del 
portico, despues que se obtiene el punto de desempeno. Se almacena en Tabla4 
el desplazamiento en el tope, el cortante basal, el coeficiente de masa modal y el 
factor de participacion modal, para cada uno de los puntos de los ciclos de carga. 
En Tabla5se almacena el primer modo de vibracion normalizado a la unidad en el 
tope en cada ciclo de carga y en Tabla6 se guarda la matriz de rigidez lateral en 
forma vectorial, solo se almacena la matriz triangular superior. 

7 ESPECTRO DE DEMANDA 

En el sistema de computacion CEINCI-LAB se tienen tres programas para 
determinar los espectros y se denominan: espectro_nec1 1_mec, que determina 
las aceleraciones espectrales de acuerdo a la Norma Ecuatoriana de la 
Construccion (2011); espectro_ern12_mec encuentra los espectros con los 
factores de sitio encontrados en el estudio de microzonificacion sfsmica de Quito 
por ERN (2012); Acel_ponderadas_mec que determina las aceleraciones 
espectrales para las cuatro micro zonas del Centra Norte de Quito, basicamente 
este programa contiene la Tabla 2 de este artfculo. 

En la tabla 9 se presentan los datos con los cuales se hallan los espectros 
indicados en la figura 21, para un lugar de Quito, ubicado en el Barrio la Floresta 
en el Centra Norte de Quito. 



Tabla 9 Datos para encontrar espectros en un sitio ubicado en la Floresta. 



Lugar 


NEC-11 


ERN-12 


ILB-14 




Suelo 




F d 


F s 


Zona 


La Floresta 


D 


1.355 


0.975 


1.225 


1 



20 




qI 1 1 1 1 1 1 

0 0 5 1 1.S 2 2.5 3 

T(s) 

Figura 21 Espectros para un lugar del barrio la Floresta que se halla en suelo tipo 
D, encontrados con NEC-11; ERN-12; ILB-14. 
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En el Metodo del Espectro de Capacidad MEC, la demanda sismica debe 
ser colocada en el formato S d - S a . Siendo S a la aceleracion espectral que en la 
figura 21 se ha denominado A y S d el desplazamiento espectral, el mismo que se 
halla a partir de la definicion del Pseudo Espectro. 



La variable no definida es T que es el periodo de vibracion. Lo que se 
desea encontrar es los desplazamientos y fuerzas laterales en la estructura de 6 
pisos de acero ante cada uno de los espectros indicados en la figura 21, aplicando 
el MEC. 



En el MEC se coloca en un grafico el Espectro de Capacidad de la 
estructura y el Espectro de Demanda; los dos espectros en el formato S d — S a . El 
MEC permite visualizar el probable desempeno de una estructura ante un sismo; 
en efecto, si los dos espectros se hallan bastante distantes la estructura va a tener 
muy mal comportamiento o va a colapsar. 

Se denomina punto de demanda, a la interseccion del espectro de 
capacidad con el espectro de demanda y que satisfagan la siguiente condicion. 
Con la demanda de ductilidad de la estructura se obtiene el parametro ficon el 
cual se halla el espectro para un factor de amortiguamiento diferente al 5%. 



Donde S dp es el desplazamiento en la interseccion de los dos espectros; 
S y es el desplazamiento cuando una seccion de la estructura ingreso al rango no 
lineal. Se considera un diagrama de histeresis bilineal, ver figura 22, para 
encontrar el factor de amortiguamiento equivalente % eq a partir de una demanda de 
ductilidad conocida fi. En Aguiar (2012) se encuentra muy bien detallado la 
determinacion de % eq la misma que fue propuesta por Jennings (1968). 




(13) 



8 METODO DEL ESPECTRO DE CAPACIDAD 



Sa 




Sap 



Sdp 



Sd 



Sdp 



Figura 22 Diagrama de histeresis bilineal. 
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E D = 2 »-!)(! -a) (15) 
eq 4 n E s n fj. (1 + afi — a) 

Donde £ D es la energfa disipada en un ciclo de histeresis; E s es la energfa 
almacenada elasticamente; a es la relacion entre la rigidez pos fluencia con 
respecto a la rigidez elastica. 

El ATC-40 al considerar las imperfecciones de las curvas de histeresis, 
en el sentido de que no son rectas como se ha considerado en el modelo bilineal 
sino curvas, introduce un factor de correccion k de acuerdo al nivel de diseho de 
la estructura. De tal manera que el amortiguamiento viscoso efectivo, es: 

W = f + kf eq (16) 

Donde <f es el factor de amortiguamiento intrfnseco de la estructura 
(normalmente 0.05). En Aguiar (2012) se encuentran los valores de k propuestos 
por el ATC-40. Ahora bien, el programa Metodo_Espectro_de_Capacidad 
determina el punto de demanda aplicando el MEC; en ese programa se 
consideran los valores de k y a indicados en la tabla 10. 

Tabla 10 Valores de k y a, que considera programa de CEINCI-LAB denominado 
Metodo_Espectro_de_Capacidad 



Nivel de Diseho 


k 


a 


L 1 


1 


0.2 


2 


0.8 


0.15 


3 


0.6 


0.10 


4 


0.333 


0.05 



El Nivel de Diseho 1, corresponde a una estructura muy bien detallada la 
misma que es capaz de tener demandas de ductilidad altas; por el contrario la que 
corresponde al Nivel de Diseho 4 tiene muy poca capacidad de disipacion de 
energfa; los otros dos niveles de diseho corresponden a valores intermedios de 
disipacion de energfa. Los valores de k son similares a los propuestos por el ATC- 
40. Se recuerda que a relaciona la rigidez post fluencia con respecto a la elastica. 

El espectro de demanda inicial es para un factor de amortiguamiento del 
5%. Ahora bien al tener un factor de amortiguamiento f e/ diferente al 5% se debe 
encontrar el parametro B que permite hallar el espectro de demanda para 
cualquier valor de amortiguamiento f e/ . 



B = 



jef_\ (17) 
0.05 



El espectro inelastico se halla dividiendo las ordenadas espectrales del 
espectro elastico para el valor de B. En el espectro inelastico es que se halla el 
punto de demanda. El calculo del punto de demanda se realiza en forma 
interactiva. 
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9 DESPLAZAMIENTOS Y FUERZAS EN CADA PISO 

El punto de demanda S dp esta asociado a un sistema de un grado de 
libertad. Para hallar el desplazamiento en el tope del edificio se utiliza la ecuacion 
(10). Pero con el factor de participacion modal FP que tiene la estructura para el 
punto de demanda; por este motivo es que se almaceno en la Tabla4 los factores 
de participacion modal. 

D t =FP* S dp 

En Tabla5 se tiene el primer modo de vibracion para cada incremento de 
carga. Ahora que se conoce en que incremento de carga se encontro S dp , se halla 
el primer modo de vibracion correspondiente que esta normalizado a la unidad en 
el tope, </>! y se multiplica por D t , con lo que se halla los desplazamientos laterales 
q en cada piso de la estructura. 

q=D t<f > 1 (18) 

Para encontrar las fuerzas laterales en cada piso, primero se encuentra el 
cortante basal V asociado con el punto de demanda, utilizando para ello la 
ecuacion (9) y despues se encuentran las fuerzas en cada piso de acuerdo al 
NEC-11. 

V = S ap ttl M T 

w x h x (19) 
F = — — — V 

x EILi^f 

T < 0.5 k = l 

0.5 < T < 2.5 k = 0.75 + 0.5 T ( 20 ) 
T > 2.5 k = 2 

Donde F x es la fuerza lateral en el piso x; w x es el peso del piso x; h x es la 
altura desde el nivel del suelo hasta el piso x; k es un coeficiente que depende del 
penodo T. 



El calculo de los desplazamientos y fuerzas en cada piso, lo realiza el 
prog ra ma despla_fuerza 
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10 RESULTADOS 

Se realizo el analisis sfsmico de la estructura de acero de 6 pisos, indicada 
en la figura 13, con cada uno de los espectros presentados en la figura 21 y los 
resultados obtenidos al aplicar el MEC, se indican en las figuras 23 a 24 para los 
espectros de: NEC-11; ERN-12; ILB-14. En todos estos graficos se presenta el 
espectro de diseno para un factor de amortiguamiento del 5% (Espectro de 
Demanda, S d - S a ) y tambien se indica el espectro inelastico con el que se obtiene 
el punto de desempeno ((5 d - S a )/B). 

En la Tabla 11 se presentan los resultados, para los tres espectros el 
periodo con el que se halla el punto de demanda es igual para los tres espectros y 
vale 1.3134 seg. Posteriormente se indica el desplazamiento y aceleracion del 
punto de demanda, S dp ,S ap . Luego se indican los desplazamientos laterales en 
cada piso desde el primero q 1 hasta el ultimo piso q 6 . A continuacion se indican las 
fuerzas laterales en cada piso y en la ultima fila el cortante basal. 
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Figura 23 Punto de desempeno con espectro de NEC-11 
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ERN12 
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0 6 Sd(m) 

Figura 24 Punto de desempeno con espectro de ERN-12 
Sa (m/s2) ILB 
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Figura 25 Punto de desempeno con espectro de ILB-14 
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Tabla 11 Desplazamientos y fuerzas laterales en cada piso encontradas con cada 



uno de los tres espectros. 



Variable 


NEC-11 


ERN-12 


ILUMBISI-LA BOTA 


T(s.) 


1.3134 


1.3134 


1.3134 


S dp (m.) 


0.1687 


0.1034 


0.1142 


m 

Sap ^s2^ 


7.0466 


5.7062 


5.9511 


Qi ( m - ) 


u.uoou 


U.UOO/ 


U.UO/ z 


q 2 (m. ) 


u. i oyo 


u.uoycs 


u.uoou 


q 3 (m. ) 


D 1704 

\J. I / U4 


n nsfin 

u.uouu 




q 4 (m.) 


0.1845 


0.1084 


0.1196 


<7s Cm.) 


0.1931 


0.1280 


0.1412 


q 6 (m.) 


0.1975 


0.1389 


0.1533 


F, (T.) 


11.2211 


7.8516 


8.6641 j 




23.5626 


16.4871 


18.1932 


F 3 (T.) 


38.2000 


26.7292 


29.4951 j 


F,(T.) 


54.6970 


38.2724 


42.2327 


Fs(T.) 


72.7877 


50.9308 


56.2010 


Pe(J.) 


52.2766 


36.5788 


40.3639 


V 0 (T.) 


252.745 


176.850 


195.15 



En la tabla 1 1 se observa que los mayores desplazamientos y fuerzas 
laterales se hallan con el espectro del NEC-1 1 , esto es debido a que en la zona de 
penodos que estan alrededor de 1.3 seg., las mayores ordenadas espectrales se 
tienen con el espectro del NEC-1 1 , ver figura 21 . 

1 1 COMENTARIOS Y CONCLUSIONES 



Varios objetivos se aspira haber cumplido con este artfculo y son los 
siguientes: presentar los resultados de una investigacion realizada en la 
Universidad de Fuerzas Armadas ESPE, para determinar espectros para el Centra 
Norte de Quito, en caso de generarse un sismo de magnitud 6.2, en el segmento 
de falla de llumbisf-La Bota ILB. 

Estos espectros presentan aceleraciones espectrales mas altas en la zona 
de penodos cortos, en un 30% aproximadamente, con relacion a los que se 
obtienen con el espectro del NEC-11 o ERN-12, pero en la zona de penodos 
intermedios y largos sus ordenadas espectrales son mas bajas (Tfpico de los 
sismos corticales) en relacion a las que se hallan con el NEC-11 o ERN-12, que a 
mas de los sismos corticales consideran la ocurrencia de sismos interface de bajo 
angulo de buzamiento que tienen varios picos caractensticos en la zona 
descendente del espectro y por ende mayores ordenadas espectrales. 

Otro objetivo que se busco con el desarrollo de este artfculo fue presentar 
en forma didactica la obtencion de la curva de capacidad sfsmica resistente que 
relaciona el cortante basal con el desplazamiento lateral maximo, mediante 
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analisis estatico no lineal; la curvas del espectro de capacidad y de demanda. 
Finalmente indicar el Metodo del Espectro de Capacidad y la determinacion de los 
desplazamientos y fuerzas laterales, en cada piso, a partir del punto de demanda. 
Toda la explicacion se complemento con el uso de los diferentes programas del 
sistema de computacion CEINCI-LAB, desarrollado en la Universidad de Fuerzas 
Armadas ESPE. 

Finalmente se realizo el analisis sfsmico de una estructura de acero de 
seis pisos con los tres espectros: NEC-11; ERN-12; ILB-14 y se obtuvo que los 
mayores desplazamientos y fuerzas laterales se hallan con el espectro del NEC-11 
debido a que el perfodo del punto de demanda se halla en la zona descendente 
del espectro donde se tienen mayores ordenadas espectrales en relacion a las 
que se producen con el espectro ILB-14. 

REFERENCIAS 

1. Abrahamson N., Silva W., Kamai R. (2013), Update of the AS08 Ground- 
Motion Prediction Equations Based on the NGA-West2 Data Set. Pacific 
Earthquake Engineering Research Center, PEER, 143 p. 

2. Aguiar R., Rivas A., Trujillo S., Gaspar J., Parra H., (2015), "Movimientos 
fuertes y espectros para el Centra Norte de Quito asociados a las fallas 
ciegas inversas", Artfculo del libro: 50 ahos de creacion de la Carrera de 
Ingenierfa Geografica, 31 p. Departamento de Ciencias de la Tierra y la 
Construccion. Universidad de Fuerzas Armadas ESPE. 

3. Aguiar Roberto, (2014), Analisis Matricial de Estructuras con CEINCI-LAB, 
Cuarta Edicion. Instituto Panamericano de Geograffa e Historia, IPGH, 
Ecuador, 676 p., Quito. 

4. Aguiar Roberto, (2013), Microzonificacion sismica de Quito, Primera 
Edicion. Instituto Panamericano de Geograffa e Historia, IPGH, Ecuador, 
212 p., Quito. 

5. Aguiar Roberto, (2012), Dinamica de Estructuras con CEINCI-LAB, 
Instituto Panamericano de Geograffa e Historia, IPGH. Segunda edicion, 
416 p., Quito. 

6. Aguiar Roberto, (2003), Analisis sfsmico por desempeho, Centra de 
Investigaciones Cientfficas. Escuela Politecnica del Ejercito, 342 p., Quito. 

7. Aguiar Roberto, (2002), Sistema de computacion CEINCI3 para evaluar 
daho sfsmico en los Pafses Bolivarianos, Centra de Investigaciones 
Cientfficas, Escuela Politecnica del Ejercito, 302 p. Quito. 

8. Alvarado A., Audin L, Nocquet M., Lagreulet S., Segovia M., Font Y., 
Lamarque G., Yepes H., Mothes P., Rolandone F., Jarrfn P., and 
Quidelleur X., (2014), "Active tectonics in Quito, Ecuador, assessed by 
geomorpholigical studies, GPS data, and crustal seismicity", Tectonics, 
AGUPLICATIONS, 17 p., Article online. 

9. Anderson, J. G. (1979). Estimating the seismicity from geological structure 
for seismic-risk studies. Bulletin of the Seismological Society of 
America, 69(1), 135-158. 

10. AISC 41 "American Society of Engineers", Seismic Rehabilitation of 
Existing Buildings, 2006. 



Roberto Aguiar, David Mora y Enrique Morales 



38 



11. ATC-40 "Seismic evaluation and retrofit of concrete buildings", Applied 
Technology Council, Redwood City, California, 1996. 

12. Beauval C, Yepes H., Bakun W., Egred J., Alvarado A., and Singaucho C, 
(2010), "Locations and magnitudes of historical earthquakes in the Sierra 
of Ecuador (1586-1996), Geophys. Journal International, 181, 1613-1633. 

13. Brune, J. N. (1968). Seismic moment, seismicity, and rate of slip along 
major fault zones. Journal of Geophysical Research, 73(2), 777-784. 

14. Campbell K., Bozorgnia Y., (2013), NGA-West2 Campbell-Bozorgnia 
ground motion model for the horizontal components of PGA, PGV and 5%- 
Damped elastic Pseudo-Acceleration response spectra for periods ranging 
fro 0.01 to 10 sec, Pacific Earthquake Engineering Research Center, 
PEER, 75 p. 

15. Chistopupoulus.C. and Filiatraul,A.(2006). Priciples of passive Suplemental 
Damping and Seismic Isolation. IUSS Press, Instituto Universitario di Studi 
Superiori di Pavia, Pavia, Pavia, IT. 

16. Chopra, A. K. y R. K. Goel (2001), "A modal pushover analysis procedure 
to estimate seismic demands for buildings: Theory and preliminary 
evaluation", Reporte No. PEER 2001-03, Pacific Earthquake Engineering 
Research Center. 

17. Chopra, A. K. y R. K. Goel (2002), "A modal pushover analysis for 
estimating seismic demands of buildings", Earthquake Engineering and 
Structural Dynamics, Vol. 31, pp. 561-582. 

18. ERN-2012, Microzonificacion sfsmica del distrito metropolitano de Quito. 
Estudio de la amenaza sfsmica a nivel local. Programa de reduccion de 
riesgos urbanos del Distrito Metropolitano de Quito. 

19. Filiatrault, A., (2014), "Structural Control", Class Notes Buffalo University. 

20. Goel, R. K. y A. K. Chopra (2004), "Evaluation of modal and FEMA 
pushover analyses: SAC Buildings", Earthquake Spectra, Vol. 20, No. 1, 
pp. 225-254. 

21. Hanks T.C. and Kanamori H. (1979), "A moment magnitude scale". Journal 
of Geophysical Research. 84: 23480-23500. 

22. Jennings J., (1968) "Equivalent viscous camping for yielding structures", 
Journal Engineering Mech. Div. ASCE, 94 ( 1 ), 103-1 16. 

23. Leonard M., (2010), "Earthquake fault scaling: Self consistent relating of 
rupture length width, average displacement, and moment release", Bulletin 
of the Seismological Society of America, 100 (SA), 1971-1988. 

24. Moroni M., editora, (2012), Mw=8.8 Terremoto en Chile, 27 de febrero de 
2010, Facultad de Ciencias Fisicas y Matematicas, Universidad de Chile. 

25. NEC-2011, Norma Ecuatoriana de la Construccion. Capitulo 2: Peligro 
Sismico y requisitos de diseho sismo resistente, 76 p. Camara de la 
Construccion de Quito. 

26. Prakash, V., G. H. Powell y Fillipou (1992), "DRAIN-2DX: base program 
user guide", Report No. UBC/SEMM-92/29, Department of Civil 
Engineering, University of California at Berkeley. 

27. Rivas A., Aguiar R., Benito M. B., Gaspar J., Parra H., (2014), 
"Determinacion del periodo de recurrencia y magnitud maxima para el 
control de las estructuras en el rango elastico ante un sismo asociado a 
las fallas inversas de Quito, Revista Internacional de Ingenierfa de 
Estructuras, 1 9 (2), 201 -21 7. 



39 



Peligrosidad sismica de Quito y el Metodo del Espectro de Capacidad 



28. Trujillo S., (2014), Espectro de Control para el Centro Norte de Quito, 
Tesis para obtener el tftulo de Ingeniero Civil. Universidad de Fuerzas 
Armadas ESPE. 

29. Wesnousky S., (1986), Earthquakes, Quaternary faults, and seismic 
hazard in California, J. Geophys, Res., 91 (B12), 12587-12631. 

30. Zhao J., Zhang J., Asano A., Ohno Y., Oouchi T., Takahashi T., Ogawa H., 
Irikura K., Thio H., Somerville P., Fukushima Y., and Fukushima Y., (2006), 
"Attenuation relations of strong ground motion in Japan using site 
classification based on predominant period", Bulletin of the Seismological 
Society of America, 96 (3), 898-913. 



